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ABSTRACT 

Aims. The X-ray emission of massive stars has been studied when these objects are in their main-sequence phase, as well as in their 
Wolf-Rayet phase. However, the X-ray properties of the transitional Luminous Blue Variable (LBV) phase remain unknown. 
Methods. Using a dedicated but limited XMM-Newton survey as well as archival XMM-Newton and Chandra observations, we 
performed the first X-ray survey of LB Vs: about half of the known LB Vs or candidate LB Vs are studied. 

Results. Apart from the well known X-ray sources eta Car and CygOB2#12, four additional LBVs are detected in this survey, 
though some doubt remains on the association with the X-ray source for two of these. For the other LBVs, upper limits on the flux 
were derived, down to log[Lx/LBOL] - 9.4 for PCyg. This variety in the strength of the X-ray emission is discussed, with particular 
emphasis on the potential influence of binarity. 

Key words. X-rays: stars - Stars: early-type - Stars : mass-loss - Stars: variables: S Doradus 



1. Introduction 

Luminous Blue Variables (LBVs), also called S-Doradus vari- 
ables, form a class of peculiar massive stars mostly found in the 
top-left part of the Her tzsprung-Russell diagram (for a review, 
see Humph reys & Davidson 1994). Their main characteristics 
are a high mass-loss rate (up to 10^"* Mq yr"'), a high luminosity 
(~ 10* Lq), and significant spectroscopic as well as photometric 
variability. These variations occur on several timescales and am- 
plitudes, from common small-scale changes of about O.lmag to 
rare giant eruptions with 2.5mag variation. 

The most emblematic objects in this category are those 
which have undergone unexpected, enormous variations (P Cyg 
in the 17th century, eta Car in the 19th century, and HD5980 
in 1993-5 and 1960-5). However, these dramatic cases should 
not be considered as the best representative objects of the class 
(iDavidsonL [19891) . 

LBVs are often surrounded by young and massive circum- 
stellar neb ulae, which result from stro ng mass-loss episodes (see 
reviews bv lLamer ill989HSmithll2011ah . 

In view of their variability and the presence of nitrogen- 
enriched circumstellar nebulae, LBVs are thought to be asso- 
ciated to a short (< 10^ yr) unstable stage in the life of massive 
stars, betw een the main-sequen ce phase and the late Wolf-Rayet 
phase (e.g. iLamers et all 1200 lb . Their crucial role in the mass- 
loss process of massive stars has been revived in the last few 
years as the mass-loss rates of O and WR stars were revised 
downwards (e.g. Smith, 2008"). 

While dedicated campaigns of X-ray observations have tar- 
geted massive stars in their O-star phase as well as in their Wolf- 
Rayet phase, not much was done for studying LB Vs. Only the X 
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ray bright sources eta Car ( Sew ard et all Il979t ICorcoran et al.L 
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I20T0I), Cy g OB2 # 12 (Harnden et all [19791: [RauwL [20TTI) and 
HD5980 jNaze et al.. 2002, 2004 l2007h . which display X-ray 
luminosities of lO"*"^ erg s"', have received attention. The sole 
exception may be P Cyg , whose ROSAT HR I data are discussed 
in length by [ Berghofe r & Wendkerl (l2000l) and mentioned by 
lOskinoval (l2005a) . This observation led to an upper limit on the 
luminosity of ~ 10^" erg s Finally, non-detections or a very 
weak detection of a few candidate LBVs are sometimes reported 
when discussing the X-ray emissi on of their h ome clusters, but 
withou t much det ail: Pistol Star (iMuno et al. I [2006; Oskinoval 
2005b'), Sher25 (IMoffat et al.L l2002h and W243 (jClark et all 
2008). 

This paper aims at bridging the gap between X-ray studies of 
O and WR stars by using the most sensitive X-ray observatories 
available at the present time, XMM-Newton and Chandra. Using 
a dedicated (though limited) survey as well as archival data, we 
investigate the X-ray emission of about half of the known LBVs. 
The paper is organized as follows: Sect. 2 presents the LBV cat- 
alog used for our survey, Sect. 3 shows the available X-ray ob- 
servations and their reduction, Sect. 4 lists our results. Sect. 5 
discusses them and Sect. 6 finally summarizes the new informa- 
tion and concludes. 

2. The sample 

Two main S-Dor or L BV c atalogs are available in the literature: 
Ivan GenderenI (l2001b and IClarketal.1 ([2005b.) . The latter lists 
more sources - 12 confirmed and 23 candidates (which will be 
noted cLBVs in the following) - and includes the 21 objects from 
van Genderen, with the exception of IRC- H 1 0420 and HP 148937 
(which is of the peculiar Of?p t ype, see IWalbornll 19721 for ini- 
tial definition of this class and N aze et al. 201 id for an in-depth 
study of this object in particular). 

Since Clark et al.'s work, two new cLBVs have been reported 
in the literature: MWC930 (Miroshnichenko et al., 2005) and 
2MASS J17460562- 2851319 (also known as GO.120-0.048, 
iMauerhan et al.l[2"010h . Both were added to our catalog. 
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Recently, Spitzer surveys of the Galactic plane (e.g. 
MIPSGAL) have led to the discovery of many ring nebu- 
lae surrounding luminou s central stars (Wachter et al., 2010; 
iGvaramadze et al.Ll2010bh and follow-up observations were con- 
ducted in order to determine the nature of these objects. 

IWachter et"an (l2010h list 4 known LB Vs (or candidates) and 
6 additional candidate LBVs. They also have 6 Be candidates 
for which they mention that their separation from the cLBVs is 
arbitrary and that these objects should therefore b e considered 
as potential LBVs too. IGvaramadze et alj (1201 Obi) reported the 
detection of 2 additional cLBVs and a (known, though not in- 
cluded in the Clark et al. compilation) cLBV. Note that one of 
these, MNl 12, was identified as Be in Wachter et al.'s list (under 
the name 2MASS J19443759-H 2419058). 

Expanding on these works, 'Wachter e t al.l (1201 1 ') list 8 ad- 
ditional LBV candidates (one having been classified as Be in 
her previous paper) and 4 additional Be candidates (two already 
classified so in their previous paper). They also mention find- 
ing 6 additional LBV candidates and one B e candidate from 
the f ollow-up observations of the sample of IGvaramadze et aP 
(I20i0b). 

Considering these new detections, we sum up Clark's list 
of 35 LBVs or cLBVs, 2 recent additions, 21 objects from 
Wachter's sample and 9 from Gvaramadze's sample. This en- 
larges the (c)LBV catalog to 67 objects. They are given in 
Table[Tl along with their main properties, when known (see 
quoted references). Typic al errors on log(r) an d log(L) are 0.05 
and 0.2 dex, respectively ("v an Genderenil200Th . More generally, 
it must be kept in mind that the properties of LBVs are much 
more uncertain than those of O and WR sta rs. Names and po- 
sitions come from the Simbad database and IGvaramadze et al.l 
(1201 Obi) for their sample - RA and DEC were rounded to 2 deci- 
mal digits. For GRS 79.29h-0.46, the Simbad and van Genderen 
coordinates do not agr ee well, and a revised estima te of the 
position was found in iJimenez-Esteban et al.l (l2010l) : the co- 
ordinates quoted in Table [T| correspond to those of 2MASS 
J20314228-H4021591, the (bright) 2MASS source closest to the 
Simbad position as well as to the van Genderen and Jimenez- 
Esteban coordinates. 



3. X-ray observations 

We have obtained XMM observations for a limited LBV sur- 
vey (ObslD=0600030, PI Naze). To enlarge the sample as much 
as possible, we further searched the archives for additional, 
serendipitous observations of the known LBVs or cLBVs lying 
within 15' (for XMM-Newton) or 10' (for Chandra) of the on- 
axis position. Only imaging exposures (i.e. no XMM "timing" 
mode or Chandra "continuous clocking" mode) longer than 5ks 
were considered. Indeed, the quality of these serendipitous ob- 
servations is often not as good as for the requested on-axis data, 
as the source usually lies far off-axis and the effective exposure 
time may not be very long. 

In total, including our targets, XMM-Newton has observed 
26 (c)LBVs while Chandra has observed 21 objects. A closer 
look indicates that 16 targets are in common, implying that 31 
(c)LBVs, or half of the catalog, have been observed in X-rays 
(see crosses in the last columns of Table [1]). Note that, if we re- 
strict ourselves to the Clark catalog, XMM-Newton has observed 
19 objects and Chandra 16, with 13 sources in common: X-ray 
observations therefore exist for two-thirds of Clark's catalog. 

Two final remarks must be made. First, the datasets for the 
two well-known X-ray emitters eta Car and Cyg OB2 #12 have 



not been re-reduced or re-analyzed, as many published papers al- 
ready discuss them in detail (see Sects. 4.1.1 and 4.1.3). Second, 
the analysis of the ten sources close to the Galactic center relies 
on the Chandra data only, as XMM-Newton is not very good at 
disentangling sources in this crowded region. 

3.1. XMM-Newton 

The XMM-Newton EPIC data were reduced in a standard 
way using SAS vlO.O. After the pipeline processing (tasks 
emproc, epproc), the recommended pattern filtering was ap- 
plied (PATTERN < 12 and XMMEA_EM filter for MOS, 
PATTERN < 4 and FLAG = for pn). The event files were fur- 
ther filtered to reject times aff'ected by background flares. Images 
in the 0.5-8.0keV domain were then produced and used for the 
detection algorithm edetectchain when only one exposure was 
availably The logarithmic likelihood for detection was put at 
10, equivalent to a Act detection. Sensitivity maps for likelihoods 
log L - 2.3 (to be consistent with the Galactic Center data, see 
Sect. 3.4) and combining all instruments were also produced us- 
ing the task esensmap. These maps provide equivalent on-axis 
limiting count rates, which were transformed into actual fluxes 
using WebPimm^. The chosen model is a simple optically-thin 
thermal emission with a temperature of 0.6 keV and an absorp- 
tion directly related to the known extinctions of the t argets (using 
A^H = 5.8 X 10^' x^rfi- V)cm ^ Bohlin et al.ll978l) . Using tem- 
peratures of 0.3 or 1.0 keV does not change the results by much: 
for an observed flux of 10"'^ ergcm"^ s"' and a representative 
absorbing column of 10^^ cm"^, the total EPIC count rates vary 
by 10% in these cases compared to the 0.6 keV case. As our lim- 
its are certainly not more accurate than 10%, it means that they 
are quite insensitive to the adopted temperature. For the Wachter 
et al. or Gvaramadze et al. sources, no information on distance, 
bolometric luminosity or reddening is available. To get a limit 
in (observed) flux, we assumed a representative absorbing col- 
umn of 10^^ cm"^. Changing the latter by an order of magnitude 
(i.e. 10^' or 10^^ cm"^) would change the total EPIC count rate 
corresponding to a given observed flux, whatever the filter, by a 
factor of ~ 2. We therefore consider that the flux limits for these 
objects are accurate within a factor of 2. 

For Chandra data, the upper limit determination is slightly 
different (see next section), but we checked on the XMM- 
Newton data of AG Car that both methods provide similar es- 
timates, and thus are fully compatible with each other 

When several exposures of a given region exist, a recent 
task added to the SAS, emosaicproc, should help combine them. 
However, its use is not straightforward. For 2MASS J 18022233- 
2238002, two of the three exposures had a strange bug in the 
attitude file which prevents their use with this task (see also be- 
low, in the discussion of this source). For W243, the start time 
of the Rev 1499 summary file was wrong and had to be cor- 
rected by hand before the task could work. In addition, while 
the detection algorithm ran otherwise smoothly, the creation of 
a sensitivity map combining all exposures and all instruments 
is not always possible, as the task esensmap only accepts com- 
binations of < 10 input files. However, the effective area is not 
always the dominant factor limiting the sensitivity of the XMM- 
Newton data (see Sects. 4.1.2 and 4.3.10) and the impact of this 
problem is thus quite limited. 



' See the analysis threads for detail on its use: 
http://xmm.esac.esa.int/sas/current/documentation/threadV] 
~ http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 
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Table 1. Catalog of the known LBVs and LBV candidates (see Sect. 2), X and C refer to XMM-Newton and Chandra observations, 
respectively. 



Simbad Name 


RA(J2000) 


DEC(J2000) 


log(reff) 


log(LBOL) 


E{B - V) 


d 


Ref 


X ? 


c? 




hh mm ss 




K 






kpc 








Clark's catalog 




















A. Confirmed LBVs 




















eta Car 


10 45 03.59 


-59 41 04.26 


4.36-4.18 


5.14-6.34 


0.50 


2.3 


vGOl 


X 


X 


PCyg 


20 17 47.20 


-1-38 01 58.55 


4.26 


5.70 


0.51 


1.7 


vGOl 


X 




HD 168607 


18 21 14.89 


-16 22 31.76 


3.97 


5.38 


1.55 


2.2 


vGOl 






AG Car 


10 56 11.58 


-60 27 12.81 


4.46-4.13 


6.14-6.22 


0.63 


6.0 


vGOl 


X 




HRCar 


10 22 53.84 


-59 37 28.38 


4.34-3.90 


5.90-5.62 


1.00 


5.2 


vGOl 






HD 160529 


17 41 59.03 


-33 30 13.71 


3.95 


5.46 


1.10 


2.5 


vGOl 


X 




WRAY 15-751 


11 08 40.06 


-60 42 51.7 


4.48 


5.91 


1.60 


6.0 


pa06,vG01 






qF362 


17 46 17.98 


-28 49 03.46 


4.05 


6.25 


8.00 


11.5 


na09 


X 


X 


AFGL 2298 


19 00 10.89 


-1-03 45 47.1 


4.01-4.07 


6.11-6.30 


9.00 


10.0 


cl09b 






GAL 024.73-1-00.69 


18 33 55.29 


-06 58 38.7 


4.08 


5.6 


3.40 


5.2 


cl03 


X 




CI* Westerlund 1 W 243 


16 47 07.50 


^5 52 29.16 


3.93 


5.86 


4.35 


4.5 


ri09,cl05 


X 


X 


GCIRS 34W 


17 45 39.73 


-29 00 26.50 


4.29 


5.5 


10.80 


7.6 


ma07 


X 


X 


B. Candidate LBVs 




















NAME VI Cyg 12 


20 32 40.96 


+41 14 29.28 


4.11 


6.42 


3.40 


1.7 


vGOl 


X 


X 


GRS 79.29-h0.46 


20 31 42.28 


+40 21 59.13 


4.40 


6.30 


5.00 


2.0 


vGOl 




X* 


Pistol Star 


17 46 15.24 


-28 50 03.58 


4.07 


6.20 


8.00 


11.5 


na09 


X 


X 


HD 168625 


18 21 19.55 


-16 22 26.06 


4.08 


5.34 


1.50 


2.2 


vG92 






zetOl Sco 


16 53 59.73 


^2 21 43.31 


4.26 


6.10 


0.66 


2.0 


cr06 






HD 326823 


17 06 53.91 


-42 36 39.74 


4.34 


5.30 


1.15 


2.0 


nir07 






HD 316285 


17 48 14.04 


-28 00 53.13 


4.18 


5.44 


1.81 


1.9 


vGOl 


X 




Hen 3-519 


10 53 59.59 


-60 26 44.31 


4.48 


6.26 


1.30 


8.0 


vGOl 






HD 80077 


09 15 54.79 


-49 58 24.58 


4.23 


6.30 


1.47 


3.0 


vGOl 






AS 314 


18 39 26.11 


-13 50 47.19 


4.01 


4.90 


0.90 


8.0 


vGOl 






MWC 314 


19 21 33.97 


+ 14 52 56.89 


4.48 


6.34 


1.84 


3.0 


vGOl 






GRS 25.5+0.2 


18 37 05.21 


-06 29 38.0 






10.00 


14.5 


vGOl 


X 




GAL 026.47-1-00.02 


18 39 32.24 


-05 44 20.5 


4.23 


6.0 


3.80 


6.5 


cl03 




X 


WRAY 17-96 


17 41 35.45 


-30 06 38.8 


4.11 


6.26 


2.84 


4.5 


eg02 






WR102ka 


17 46 18.12 


-29 01 36.6 


4.40 


6.5 


8.0 


8.0 


ba08 


X 


X 


NAME LBV 1806-20 


18 08 40.31 


-20 24 41.1 




6.3 


9.4 


11.8 


fi04,ei04 


X 


X 


NAME Sher 25 star 


11 15 07.8 


-61 15 17 


4.35 


5.9 


1.60 


6.3 


sm02 




X 


[OMN2000] LSI 


19 23 47.64 


+14 36 38.4 


4.12-4.14 


5.75 


3.50 


6.0 


cl09a 




X 


GCIRS 16NE 


17 45 40.26 


-29 00 27.09 


4.29 


5.9 


8.30 


7.6 


ma07 


X 


X 


GCIRS 16C 


17 45 40.13 


-29 00 27.64 


4.24 


5.9 


8.30 


7.6 


ma07 


X 


X 


GCIRS 16SW 


17 45 40.12 


-29 00 29.06 












X 


X 


GCIRS 16NW 


17 45 40.05 


-29 00 26.87 












X 


X 


GCIRS 33SE 


17 45 40.02 


-29 00 31.0 


4.26 


5.75 


9.00 


7.6 


ma07 


X 


X 


New additions (cLBVs) 




















2MASS J17460562-2851319 


17 46 05.625 


-28 51 31.92 




6.2-6.6 


8.00 


11.5 


mulO 


X 


X 


MWC 930 


18 26 25.24 


-07 13 17.8 


4.34 


5.5 


2.5 


3.5 


miOS 


X 




Wachter's sources (cLBVs) 




















2MASS J15484207-5507422 


15 48 42.07 


-55 07 42.21 
















2MASS 116290377-4746264 


16 29 03.78 


-A7 46 26.48 
















2MASS J 16364278^656207 


16 36 42.78 


-46 56 20.73 














X* 


2MASS J17082913-3925076 


17 08 29.14 


-39 25 07.68 
















2MASS J171 10094-3945174 


17 11 00.94 


-39 45 17.45 












X 




2MASS J18415965-0515409 


18 41 59.65 


-05 15 40.93 














X 


2MASS J16431636-4600424 


16 43 16.37 


^6 0042.42 
















2MASS J16493770-4535592 


16 49 37.70 


^5 35 59.27 
















2MASS J17435981-3028384 


17 43 59.85 


-30 28 38.5 
















2MASS J18133121-1856431 


18 13 31.21 


-18 56 43.20 
















2MASS J19325284+ 1742303 


19 32 52.85 


+ 17 42 30.33 
















2MASS J19443759+2419058" 


19 44 37.60 


+24 19 05.87 










gvlO 






2MASS 116313781^814553 


16 31 37.82 


^8 14 55.30 
















2MASS J16461734^508478 


16 46 17.35 


^5 08 47.85 












X 




2MASS J17374754-3137333'' 


17 37 47.54 


-31 37 33.38 
















2MASS J17374730-3137370'' 


17 37 47.31 


-31 37 37.07 
















2MASS J17421401-2955360 


17 42 14.02 


-29 55 36.06 
















2MASS J18022233-2238002 


18 02 22.34 


-22 38 00.24 












X 




2MASS J18455593-0308297 


18 45 55.94 


-03 08 29.72 
















2MASS J18510295-0058242 


18 51 02.95 


-00 58 24.21 
















2MASS J 180705 16-20 15 163 


18 07 05.17 


-20 15 16.31 
















2MASS J19011669-H0355108 


19 01 16.69 


+03 55 10.81 
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Simbad Name 


RA(J2000) 


DEC(J2000) 


log(reff) 


log(LBOL) 


E(B -V) d Ref 


X ? 


C ? 




hh mm ss 




K 


Le 


kpc 






Gvaramadze 's 


sources (cLBVs) 














MN13 


13 42 33.08 


-62 48 11.3 












MN jy 


lO lU /0.35 


CI ^ 1 TCI 

-51 zl 25.3 












MN41 


16 26 34.28 


-50 21 01.9 












MN53 


17 09 24.78 


-40 08 45.6 








X 


X 


MN79 


18 28 33.41 


-11 46 44.2 








X 


X 


MNlOl 


19 06 24.54 


08 22 01.6 












MN107 


19 24 03.34 


13 39 49.4 












MN46 


16 43 16.37 


-46 00 42.4 












MN83 


18 39 23.01 


-05 53 19.9 










X 



Reference s are: b a08 = iBamiske et a"D 120081), cl03 = IClark et al.' (I 20q3h, cl 05 = |Clark et al.' ('2005al), cl09a HClark et al.l (12009 j), cl09 b 
= IClarketalJ (l2009bh. cr06 = Crowther et al. (200 ^), eg02 = Egan et alT T2002h. ei04 = Eikenberrv et al. ( 2004 ), fi04 = iFiger et alT d2004h . 



gvlO = iGvaramadze et al.' /2010ai). ma07 = J Vlartins et alj J2b07i). mr07 = | M arcolino et al. (2007 ), mulO = iMauerhan et al .' ^20l'^mi05 = 
^iroshnichenko et al. (2005), na09 = lNaiarro et al. (2009), pa06 4Pasauali et al. ( 2006) , ri09 = Ritchie et aTI ( l2009l) . sm02 = ISmartt et al. ( 200j) . 
lyan Genderen et alj 11993) , vGOl = lvan Genderei] J2OOII) . 



vG92 

Not yet publicly available 
'"^^ In Gvaramadze et al. ( 2010b), this source appears as MNl 12. 

*'^*IWachter et al. (201 1) quote as LBV candidate their source #27 of IWachter et alj ( l2010l) . However, in the latter paper, there are two #27 sources, 
a and b. Both are thus quoted here together. 



Table |2] lists the results of our analysis of the XMM-Newton 
data. 

Finally, a check was made using the XMM-Newton upper 
limit servefl In 3 cases, only slew data were used, leading to 
very large upper limits compared to data from pointed observa- 
tions, as could be expected. In addition, several detections with 
varying count rates were reported for W243 as well as a detec- 
tion for P Cyg, but they are due to source confusion in these re- 
gions (see details on these sources below). In all cases apart from 
those just discussed, the server finds upper limits which are com- 
patible with ours, though slightly higher. 

3.2. Chandra 

The pipeline data (level 2) available in the Chandra archives 
were downloaded but not re-processed. Further work was done 
with CIAO 4.2 and CALDB 4.3. L We first defined the source 
and background regions. The sources regions are circular, with 
radii of 5" (i.e. lOpx), while the background regions are chosen 
as surrounding annuli of radii 5" and 15". The sources W243 and 
LBV 1806-20 have close X-ray companions, so that the annular 
background regions were replaced by nearby circular regions of 
5" radii in these two cases. We then used the task specextract to 
get unbinned spectra of source and background regions, as well 
as their corresponding response matrices. When several observa- 
tions are available, the task combine spectra was used to co-add 
the spectra. 

If the source is detected, the flux in the 0.5-8.0keV band 
is then readily estimated using an absorbed optically-thin ther- 
mal model, with temperature fixed to 0.6 keV whenever less than 
lOOcts are recorded and free otherwise, and absorptions being 
always fixed to their known values (see previous section). In case 
of non-detection, the count numbers in the two regions, their 
areas, and the exposure times provide upper limits with 90% 
significance using the task aprate^ which relies on Bayesian 



' http://xmm.esac.esa.int/UpperLimitsServer/ 

* See http://cxc.harvard.edu/ciao/threads/upperlimit/' - we used it 
with alpha=0.99 (the PSF fraction encircled by the source region) and 
beta=0.01 (the PSF fraction encircled by the background region). 



statistics to evaluate errors. Using the calculated response matri- 
ces and a "fakeit" run under Xspec, we can then transform these 
count rates into actual fluxes assuming the same model as above. 

In two cases, the available data were taken using HRC rather 
than ACIS, implying that no spectral information is available. 
For LBV 1806-20, two other ACIS datasets are available. We 
thus only checked that the limit derived in the same regions using 
the task aprates on the HRC data was compatible with the ACIS 
results. For MN53, the HRC data are the sole available, and the 
HRC limit on the count rate found with aprates was transformed 
into actual flux using WebPIMMs and the above model. 

Table[3]lists our results. As for XMM-Newton, we consider 
to have a secure detection when the source is detected with at 
least 4 cr significance. 

Finally, a search was made for the cataloged (c)LBVs in the 
Chandra Source Catalog (CSCQ X-ray sources were only found 
close to GAL 026.47 -hOO.02, W243 and Sher 25. For the former, 
the X-ray source lies at 0.2", implying detection, but for the last 
two objects, the CSC X-ray sources lie at > 5", casting doubt on 
their identification with the LB Vs. 

3.3. Proprietary data 

Sher 25 is a special case in our survey since 5 observations ex- 
ist but 4 of these were still proprietary at the time of the anal- 
ysis. These additional data greatly enhance the detection limit, 
though, as the total exposure time is then multiplied by 10. The 
PI of these data, L. Townsley, kindly made them available to 
us. The full dataset (the public exposure and the 4 private ob- 
servations) was reduced and combined using the ACIS Extract 
software packag^by P. Broos: we use here an extraction specif- 
ically made at the position of Sher 25 on the full dataset. 

Similarly, a deep analysis of the X-ray emission of MCW 
930 was possible because the PI of the XMM-Newton data, L. 
Bassani, kindly made the observation available to us before the 



http ://cxc . harvard. edu/c sc/index. html 

The ACIS Extract software package and User's Guide are avail- 
able at http://www. astro.psu.edu/xray/acis/acis-anaiysisiitiiiil , see also 
iBroos etalj(l2010h . 
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Table 2. XMM-Newton imaging observations of (c)LBVs. A null off-axis angle (OFA) indicates an on-axis observation. 



Simbad Name 


Obsid 


Rev 


OFA 

(') 


EPIC count rate 

ct s"' 


erg cm ' s" 




erg S"' 


log[Lx/LBOL] 


PCyg 


0600030201 


1891 


0. 


<0.0010 


< 1.0 X 10" 


13 


< 8.3 X lO-** 


< -9.4 


AG Car 


0600030101 


1853 


0. 


<0.0009 


< 8.1 X 10" 


16 


< 1.0 X 10^' 


< -8.7.. 8.8 


HD 160529 


0600030701 


1797 


0. 


<0.0008 


< 7.8 X 10" 


16 


< 3.2 X 10'° 


< -8.5 


GAL 024.73+00.69 


0301880301 


1145 


5.6 


<0.006 


< 8.4 X 10" 


15 


< 7.2 X 10^2 


< -6.3 


CI* Westerlund 1 W 243 


0410580601 


1317 


0.5 


<0.0017 


< 2.6 X 10" 


15 


< 2.3 X 10^2 


< -7.1 




0505290201 


1409 


0.5 














0505290301 


1499 


0.5 














0555350101 


1593 


0.5 














0604380101 


1778 


0.5 












HD 316285 


0112970101 


0145 


7.8 


<0.0013 


< 1.4 X 10" 


15 


< 6.1 X 10'" 


< -8.2 


GRS 25.5+0.2 


0400910301 


1256 


10. 


<0.0025 


< 5.7 X 10" 


15 


< 1.7x lO''' 




NAME LBV 1806-20" 


0205350101 


0869 


0. 


<0.01 


< 2.2 X 10" 


14 


<4.1 X lO'"* 


< -5.3 




0164561101 


0884 


0. 














0164561401 


1066 


0. 














0406600301 


1157 


0. 














0406600401 


1237 


0. 














0502170301 


1428 


0. 














0502170401 


1523 


0. 














0554600301 


1601 


0. 














0554600401 


1691 


0. 












MCW 930* 


0650591501 


2059 


7.5 


<0.0006 


< 1.9 X 10" 


15 


< 5.0 X 10" 


< -7.4 


2MASS J171 10094-3945174 


0502080301 


1431 


0.9 


<0.003 


< 3.2 X 10" 


15 






2MASS J16461734-4508478 


0164561001 


0873 


2.6 


<0.004 


< 4.9 X 10" 


15 






2MASS J18022233-2238002 


0135742801 


0600 


11. 


<0.0035 


< 3.8 X 10" 


15 






MN53* 


0148690101 


0681 


6.9 


<0.0005 


< 1.5 X 10" 


15 






MN79 


0051940401 


0229 


3.7 


<0.0012 


< 1.3 X 10" 


15 






eta Car 


see Corcoran et ah (2010) 


0.06..3 X 10 


-10 




—5 



NAME VI Cyg 12 see Rauw (201 D 2.7 x 10"'^ 8.2 x 10" -6.1 



Notes. Column 1 provides the (c)LBV name, columns 2 to 4 the details of the observation(s) used, and the next columns list the derived EPIC 
count rates, fluxes, and luminosities in the 0.5-8.0 keV energy range. Limits correspond to a 90% chance of being detected if brighter. 

Three additional observations (Obsid - Rev = 0148210101 - 0607, 0148210401 - 0701, and 0164561301 - 0960) are available but are totally 
affected by strong flares: they were thus discarded from the analysis. Refers only to EPIC MOS data. 



end of her proprietary period. These data were reduced as our 
other XMM-Newton data (see Sect. 3.1). Note that there is also 
a Chandra dataset (400638 - 7528) for this source, but it has an 
exposure time of only 5ks ; moreover, this cLBV appears on the 
edge of a CCD, prohibiting the derivation of any meaningful flux 
limit. 

It may be worth noting that 2 additional sources (GRS 
79.29+0.46 and 2MASS J 16364278-4656207, see Table© have 
Chandra data, but they are not yet publicly available. 

3.4. The Galactic Center survey 

Ten of the cataloged cLBVs lie close to the Galactic Center A 
su rvey of this region, using 2Ms of Chandra time, was reported 
bv lMuno et alJ (l2009l) . This paper and its associated on-line re- 
source£]provide source lists and source spectra, as well as sensi- 
tivity maps in the 0.5-8.0 keV energy range (see their Fig. 9). We 
use the sensitivity map where sources brighter than the quoted 
limiting photon flux have > 90% chance of being detected (the 
other one available corresponds to a > 50% chance). This 90% 
probability corresponds to 1 .65cr or a logarithmic likelihood of 
2.3. 



' See http://www.srl.caltech.edu/gC-project/xray.h tniTI and tables on 
the journal webpage (CDS tables are incomplete and tables on the 
Caltech website contain errors). 



4. Results 

In this section, we present the results of this survey. We begin 
by presenting the detections, and then continue with the doubt- 
ful detections and the non-detections. Whenever the X-ray emis- 
sion of the object has been discussed in the past, it is clearly 
mentioned. 

4.1. Detections 
4.1.1. eta Car 

Of all LB Vs, eta Car is certainly the most spectacular, luminous, 
and best studied object... but not nec essarily the best understood 
one (Davidson & Humphreys^ Il997h ! Its fame comes from the 
two impressive brightenings observed in the 19th century (the 
"Great Eruption" in mid- 19th century and the smaller event of 
1890). Since then, there has been a lot of discussion about the 
variability of this star, and its potential recurrence times cale. A 
period of 5.5 yr was finally identified bv lDamineUl (Il996l) . and it 
is now detected from radio to X-ray energies. 

Though there is no direct evidence for a companion, this 
timescale is often interpreted as the orbital period of an ec- 
centric binary (for a review see Corcoran, 2011). Indeed, a 
colliding-wind binary (CWB) can easily explain many properties 
of the observed variability (e.g. the nebula's changes are related 
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Table 3. Chandra imaging observations of (c)LBVs. A 'H' after the ObsJD indicates HRC data, the off-axis angle (OFA) is the 
distance to the non-offset aimpoint. 



Simbad Name 


seqnum 


ObsJD 


OFA 


count rate 






l0g[Lx/LB0L] 








(') 


ct s 


-2 -1 

erg cm s 


erg s"' 




CI" West. 1 W 243 


200344 


5411 


1.7 


0.00029. .61±0.00015: 


1.7. .3. 5 X 10 : 


1.5. .3.1 X 10 : 


-7.3.. - 7.0: 




200344 


6283 


2.1 










GAL 026.47+00.02 


400603 


7493 


4.3 


0.0161+0.0009 


(1.73 ± 0.10) X 10 


(4.3 ± 0.2) X 10 


-5.95 ± 0.02 


NAME LBV 1806-20" 


500599 


625 IH 


0.4 


<0.0009 


< 3.7 X lO"''* 


< 6.9 X 10'" 


< -5.0 




500042 


746 


0.6 


<0.0011 


< 1.0 X 10"'"* 


< 2.0 X 10''" 


< -5.6 


NAME Sher 25 star 


200058 


633 


0.6 


<0.000017 


< 1.1 X 10 


< 4.9 X 10'" 


< -8.8 




ZUUODO 


LZjZo 


U.o 












ZUUODO 


12329 


0.5 














12330 


0.5 












ZUUDDD 


13162 


0.5 










rriA/TM9nnm t q i 

[UlVllN ZVJUU J 1 


9001 ^9 


2524 


5.9 


<0.00018 


^ 1 zL V 1 


^ 1 A V 1 0^2 


^ 7 1 




9001 ^9 


3711 


5.9 














7552 


6.2 


<0.00016 


< 1.3 X 10"" 






-05 1 540Q 
















MN53 


500140 

^ \J\J ± 


1963H 


7.4 


<0.002 


< 2.4 X 10~'* 






MN79 


500629 


6675 


1.7 


<0.00025 


< 1.0 X 10 






MN83'' 


400603 


7493 


5.2 


<0.0002 


< 1.5 X 10"'^ 






CLa v-dl 




see Corcoran et al. (2010) 


n Of, T, s( 10"'" 




J 


NAMF VT Cvo 1 2 






see Rauw (2011) 


2.7 X 10"'- 


8.2 X 10'' 


-6.1 


Cjcildctic Center 








photon fluxes (ph cm '^s ') 








qF 362 








< 1.5 X 10 


< 5.6 X 10"'^ 


< 7.8 X 10" 


< -5.9 


GCIRS 34W 








(9.8 ± 0.2) X 10~* 


7.5. .8.1 X 10"''' 


1.0.. 1.1 X 10" 


~ -6.1 


Pistol Star 








< 1.5 X 10 « 


< 5.6 X 10"'^ 


< 7.8 X 10" 


< -5.9 


WR102ka 








< 2 X 10^* 


< 7.4 X 10"'^ 


< 5.0 X 10" 


< -6.4 


GCIRS 16NE 








< 3.7 X lO^*^ 


< 1.4 X 10"'" 






GCIRS 16C 








< 3.7 X 10^'' 


< 1.4 X lO"'-* 


< 9.1 X 10" 


< -5.5 


GCIRS 16SW 








< 3.7 X lO^*^ 


< 1.4 X 10"'" 






GCIRS 16NW 








< 3.7 X lO^*^ 


< 1.4 X 10"'" 


< 9.1 X 10" 


< -5.5 


GCIRS 33SE 








(1.41 ±0.02) X 10"^ 


(1.42 ± 0.02) X 10"'^ 


(1.23 + 0.02) X 10" 


-6.244 ± 0.006 


2MASS J17460562 








< 1.2 X lO^*^ 


< 4.4 X 10"'^ 


< 6.2 X 10" 


< -6.0..6.4 


-2851319 

















Notes. Column 1 provides the (c)LBV name, columns 2 to 4 the details of the observation used, and the next columns present the derived count 
rates, observed fluxes, and absorption-corrected luminosities in the 0.5-8.0 keV energy range. For the detected objects, the quoted errors correspond 
to 1-cr; for non-detections, the quoted limits correspond to a 90% chance of being detected were the source brighter. 

*"* One additional observation (seqnum=500598, obsid=6224) covered the region containing this cLBV, but at a very large off-axis angle: the 
presence of a neighbouring bright pulsar prevents any detection at this position. 

Two additional observations (seqnum=500759, obsid=7630 and 9754) covered the region containing this cLBV, but the source lies on the CCDs 
outer edges. 

sion. The X-ray Ughtcurve is however not strictly periodic. The 
2003.5 event was more luminous than the 1998 event, and the 
recent 2009 event displayed a lower and softe r flux before pe- 
riastron, as we l l as a shorter minimum state ( CorcoranL 120051 : 
ICorcoran et a H I20T0I) . This may be linked to changes in the 
wind density, thro ugh variations o f the wind velocity, mass-loss 
rate or both ( Corcoran et al.L 12010) . The observed luminosities 
(in the 2-10. keV range) vary from 2 x 10^^ erg s~^ at maximum 
to value s 50 to 100 times smaller at minimum jlshibashi et all 
1999; Co rcoran et al.. 2010). As the actual value of the absorb- 
ing column and bolometric luminosity are still debated, it is dif- 
ficult to give accurate log[Lx/i'BOL], but eta Car may certainly 
reach values of -5. An example of its X-ray spectrum is shown 
in Fig.[T] 



to photo-excitation variations linked to the collision effects), 
though an exact modelling reproducing all obser ved details is 
gener ally difficult (see e.g. in X-rays the work of iParkin et al.L 
I2011h . 

In the X-ray domain, it is quite remarkable that eta Car has 
been one of the firs t massive stars detected at these energies 
(ISewardetallll979h . soon after a few Cygnus OB2 stars (see 
Sect. 4.1.3 below). It has been intensively studied since then, 
and the observed hardness of the emission and its variabil- 
ity properties have indicated a strong similarity with the well- 
known CWB WR 140 (for recent campaigns on these objects, 
see ICorc oran. 201 1; Williams, 201 1). A n exten sive RXTE mon- 
itoring fcorcoran, ,2005: Corcoran et all 1201 Oh notably showed 
that the X-ray emission appears rather stable around apastron, 
but rises before periastron, as the two stars come closer together 
At periastron, an abrupt decrease occurs, and this short mini- 
mum state is often interp reted as an occulta tion or eclipse (com- 
bined with a switch-off. iParkin et ani201 II) of the CWB emis- 
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4.1.2. W243 



Using two Chandra observations. IClark et alj(l2008l) studied the 
X-ray emission of stars in the Westerlund 1 cluster These au- 
thors mention a "weak detection" for W243: their Table 4 lists 
7 + 6 net counts for this object, corresponding to a luminosity 
of about 10^*^ erg s"' according to their count-to-flux conversion 
rat^l 

We re-analyzed the same Chandra data, with a different 
technique (see Sect. 3.2) and we also take a look at the many 
XMM-Newton datasets encompassing the source. 

As the Chandra data show that there is a source at 7.5" south 
of W243 (CXOU J164707.6-455235), we do not use an annu- 
lus to estimate the background but rather a nearby circle devoid 
of sources (see Sect. 3.2). When extracting the data at the po- 
sition of W243 in both Chandra datasets and combining them, 
we found 34 net counts for the source if the source region has 
5" radius, and 15 net counts if we reduce the extraction radius 
to 3.4" (to avoid contamination by the neighbour). This leads to 
count rates of 0.0006 1 ± 0.000 1 5 and 0.00029 + 0.000 1 cts s" ' , 
respectively. Both results are compatible with each other (e.g. 
count rates are within 3cr), but only the first case constitutes 
a formal detection with our definition (see Sect. 3.1 and Fig. 
[1]). The region used in the first case is also the one most prone 
to contamination from the neighbouring object. Indeed, chang- 
ing the extraction radius to 3.4" yields similar results for other 
Chandra targets. We therefore consider the count rate as uncer- 
tain, advocating for a "weak detection" status at best. 

All XMM-Newton data were taken with the medium fil- 
ter except for the first two datasets (0404340101 - 1240 and 
0311792001 - 1243) which we therefore exclude from the mo- 
saicking process. There is no trace of the source due south of 
W243. In fact, this source and W243 both lie on the PSF wing 
of the pulsar PSR 11647^552 which is situated at 31" east of 
W243. The presence of this bright X-ray source severely limits 
the sensitivity achieved at the position of W243. Therefore, since 
the LBV seems undetected in the XMM-A^ewto« data, an upper 
limit on its flux was derived by adding a fake source, with a sim- 
ilar PSF, at the position of the LBV|_|. We found that an X-ray 
source ~ 400 times fainter than the pulsar's emission would be 
barely discernible, implying an upper limit on the EPIC count 
rate of 0.0017 ctss"' (Table |2]i. Note that, as the pulsar was 
brighter in the two excluded observations, we obtained a better 
sensitivity limit here than when considering the whole XMM- 
Newton set. The limit derived from XMM-Newton data is quite 
similar to the values found using the weak Chandra detection, 
which means that XMM-Newton is on the verge of detecting 
W243. However, the Chandra observatory appears better suited 
for the purpose of studying W243 - deeper observations would 
help resolve the uncertainties on its properties. 

In any case, the modest X-ray emission level of 10-^^ erg s 



obser v ed in W243 is rather typical of OB stars dHarnden et alj 
1 1 979 1 'Ku & Chanan' '19791 see also recent cluster studies 
ofHana et al. 2006; AntokhineLalJ l2()08t iNaze etalJl2011bh 
as well as its l og[Lx/LBOLl of -7 ( Pallavicini et all 119811; 
iBerghoefer et all, 119971; iNazel l2009l l2011ah . This leads to two 
possibilities. On the one hand, W243 could be a single LBV with 
a wind comparable to that of a "normal" O-type star While its 



see their Table 4; this probably corresponds to a conversion towards 
absorption-corrected flux, though it is not clearly indicated in the paper. 

' The astrometric errors of the XMM-Newton datasets are negligible 
as the error on the pulsar's position, found by comparing the Simbad 
coordinates with those in the detection list of the combined mosaic 
amounts to only 0.2". 



mass-loss rate is indeed similar to that of O-type stars, the wind 
terminal velocity of W243 was found to be only 165 k m s~', 
i.e. one order of magnitude lower than in O-stars ( Ritchie et al.L 
2009). Such a slow wind casts doubt on the "LBV wind-shock 
model" scenario, as the (potential) intrinsic shocks will be rather 
weak to generate enough X-rays, compared to O-stars. On the 
other hand, W243 could actually be a multiple system composed 
of an LBV and an O-type star whose intrinsic high-energy emis- 
sion is the actual source of the recorded X-rays, rather than the 
LBV itself. 



4.1.3. NAME VI Cyg 12 = Cyg OB2#12 

With Mv ^ -10, Cyg OB2 #12 is consider ed to be one of the op- 
tically brightest stars in our Galaxy (e.g. iMassev & ThompsonI 
1991). Despite this high luminosity placing it above the 
Humphreys-Davidson limit, Cyg OB2 #12 lacks some of the 
LBV characteristics, such as temperature changes, enrichment, 
or traces of circumstellar material (Clark et al., 201 1). This ob- 
ject was one of the very first massive stars detected as an 
X-ray source during an EINSTEIN observation of CygX-3 
(iHarnden et al.L[T979l) . This star has subsequently been observed 
by all major X-ray satellites. Most recently, R auwl ( 1201 1 ) per- 
formed an in-depth monitoring of this object using six XMM- 
Newton pointings obtained in October-November 2004 and 
April-May 2007. The EPIC spectra of this star (Fig . [TJ ap- 
peared quite hard: it was fitted with an absorbed two-temperature 
thermal plasma model consisting of plasma at temperatures of 
kT ^ 0.75 and 2keV (where the strength of the second compo- 
nent is quite similar to that of the first one, not much smaller). 
The average absorption-corrected X-ray flux in the 0.5-8.0keV 
domain was 2.4 x 10"" ergcm"^ s"', which corresponds to an X- 
ray luminosity of 8.2xl0''-' erg s"' assuming a distance of 1.7kpc 
and log[Lx/i'BOL] = -6.1 (using reddening, distance an d bolo- 
metric luminosity of Table[T]rather than the values used bv lRauwl 
1201 II - though the final log[Lx/i'BOL] value are then equal). 

These temperatures, especially the second one, cannot be ex- 
plained by shocks in a slow (150km s~' ) moving wind as pro- 
posed by Kolchko va & Chentsovl (l2004l) and more rece ntly by 
.Clark et all (1201 1, 400km s ^). However. ISoiiza & Lutd (ll98oi) 
reported a variable absorption bluewards of Ha that they tenta- 
tively interpreted as a blueshifted Ha absor ption arising in an ex- 
panding shell at a velocity of 1400 km s ' . Leitherer et al.l(ll982l) 
actually inferred this velocity to correspond to the wind veloc- 
ity of CygOB2 #12. An expanding shell of even higher veloc- 
ity (around 3100 kms"') was reported by Wang & Zhu (2003^ 
A faster wind than inferred bv lKolchkova & Chentsov (2004) is 
therefore possible, though the tremendous interstellar reddening 
unfortunately hampers a reliable determination of the star's wind 
velocity. 

On the other hand, the EPIC data further revealed X-ray flux 
variability (with an amplitude of ~ 10%) on timescales rang- 
ing from a few days to a few weeks, as well as larger variations 
(37%) on timescales of years (iRauwl |201 ih . The combination 
of these properties (strong, variable X-ray emission, including a 
strong high-temperature component) suggest that Cyg 0B2 #12 
could be either an eccentric colliding wind binary (with a sec- 
ondary star featuring a fast wind) or a single star featuring a 
strong magnetic field that confines the stellar wind near the mag- 
netic equator 
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0.5 1 2 5 10 

channel energy (keV) 

CI* Westerlund 1 W 243 




12 5 
Eneroy (keV) 



Cyg OB2 #1 2 (3rd XMM obs) 




12 5 10 

Energy (keV) 




Fig. 1. X-ray spectra of the detected LBVs (from top to bot- 
tom: et a Car - frorn the 2XMM database, W243, Cyg OB2 #12 
- from ' Rauwl 12017 . and GAL 026.47+00.02), along with their 
best-fit model (see text for detail). The eta Car spectrum corre- 
sponds to XMM-MOS2 data taken at the June 2003 maximum, 
which are dominated by the soft emission of the outer ejecta 
below 1.5 keV and by c olliding-wind emisson above that limit 
dHama^uchi et al.L 120071) . The W243 spectrum corresponds to 
the 5" extraction region, and was re-binned by a factor of 10. 
^ote the strong iron line at 6.7 keV for eta Car and Cyg OB2 



4.1.4. GAL 026.47+00.02 

GAL 026.47+00.02 was detected by Chandra as a rather faint 
X-ray source (only ~ 300 counts recorded in 20 ks). This X-ray 
source is well isolated, thereby giving credit to its association 
with the cLBV. 

A rough (in view of the S/N) X-ray spectrum could be ex- 
tracted for this object (Fig.[TJ. A fit by an absorbed optically thin 
thermal plasma model yields a temperature of 0.47 + 0.11 keV 
and a normalization factor of (5.7±0.4)x lO^^cm"^ when the ab- 
sorption is fixed to the ISM value. The resulting log[Lx/i'BOL] = 
-6 is similar to that of Cyg OB2 #12, though the former object 
appears somewhat softer Such a hard and bright X-ray emission 
suggests that there are colliding winds or magnetically confined 
winds in the system. However, follow-up observations (e.g. to 
test whether the X-ray source varies, and how) would be required 
to further consti'ain the nature of GAL 026.47+00.02. 



4.2. Doubtful detections 
4.2.1. GCIRS 34W 

This LBV appeai-s at 0.3" of the X-ray source CX OUGC 
J1745 39. 7-290026 in the Galactic Center survey ( Munoetali 
I2OO9I) . However, the X-ray source does not appear point-like but 
rather looks like a diffuse emission region. This casts doubt on 
the identification of the source of the recorded X-rays with the 
sole LBV. 

This conclusion is supported by the appearance of the X-ray 
spectrum (Fig. |2]i. An absorbed optically-thin thermal plasma 
provides a rather good fit but only with extremely high temper- 
atures: kT > 53keV, norm = (1.09 + 0.05) x 10"'*cm-\ and 
an absorption fixed to the ISM value yield - 1-5. An ab- 
sorbed powerlaw is clearly better: F = 0.95 + 0.16, norm = 
(1.1 ±0.2) X lO"*" cm"^, and an absorption fixed to the ISM value 
yield = 1.1. Both models yield quite similar flux values, 
though (both are given in Table |3]l. 



4.2.2. GCIRS 33SE 

The cLBV is at 0.3" of the X-ray sour ce CXOUGC J174540.0- 
290030 of the Galactic Center survey dMuno et al.Lr2009 ). This 
source was found to vary in flux by two orders of magni- 
tude (from photon fluxes of 1.6 x lO^^phcm"^ s"' to 2.3 x 
lO-^hcm-^s-i). 

The mean spectrum appears extremely hard , with all pho- 
tons at energies larger than 2keV (Fig. |2]i. It could not be fit- 
ted by an absorbed powerlaw, an absorbed optically-thin thermal 
plasma, an absorbed bremsstrahlung, or an absorbed disk black- 
body model. In fact, a good fit was only achieved by adding a flat 
powerlaw to a line (either a gaussian or a relativistic one): for the 
lattercase (model diskline in Xspec), we found F - -0.08+0.12, 
normpL = (2.2+0.6)xl0-''cm-^ = 6.41+0.05keV, / > 77°, 
and normiine - (4.8 + 0.4) x 10"'' cm"^ 

This type of spectrum and variations are rather atypical of 
massive stars (even considering the "hard" colliding-wind sys- 
tems) and would better fit in the context of X-ray binaries. 
Therefore, we doubt that the survey has detected the true emis- 
sion of GCIRS 33SE. 
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Fig. 2. X-ray spectra of the X-ray sources apparently associated 
with GCIRS 34W and GCIRS 33SE, along with their best-fit 
model (see text for detail). Note the absence of lines, especially 
in comparison with Fig.[T] 

4.3. Non-detections 
4.3.1. P Cyg 

The ROSAT data of P Cyg wer e ana lyzed by 
iBerghofer & Wendked (l200d) and lOskinoval (l2005a). However, 
these X-ray data suffer from strong contamination by the 
intense UV emission from this hot star, which implies most 
probably a non-detection of PCyg in this dataset. Limiting 
(absorption-corrected) luminosities of 9 x 10""' erg s"' and 
2 X 10^' erg s'' in the ROSAT range (0.1 -2.5keV) were re - 
ported by IBerghofer & Wendked (l2000h and lOskinoval (l2005al) . 
respectively. This corresponds to log[Lx/iBOL] below ~ -8 
(considering the same soft energy range). 

The new XMM-Newton observation was taken with the thick 
filter, to avoid any UV contamination. In the new data, three faint 
X-ray sources are clearly detected towards the east of PCyg. 
They are at distances of 28" (for two of them) and 62" (for the 
last one of the trio). This is well above the potential astrometric 
errors and the PSF size, i.e. these source are too distant for being 
reasonably associated with the LBV. 

Definitely, there is no X-ray source at the position of P Cyg. 
The upper limit on the flux is more than an order of magnitude 
lower than the ROSAT limit though it corresponds to a harder 
energy range (0.5-8. keV for XMM-A^ewfo« vs 0. 1-2.5 keV for 



ROSAT). The signal recorded by ROSAT was thus clearly dom- 
inated by UV contamination. 

4.3.2. AG Car 

The field of AG Car looks very similar to that of P Cyg: only faint 
X-ray sources scattered everywhere, without anything at the po- 
sition of the LBV. Here, the closest X-ray source lies at nearly 
1 ' of AG Car, and there is thus no doubt on the non-detection of 
the LBV. 

Note that Hen 3-519 is rather close to AG Car, but no con- 
straint on its X-ray flux can be derived from these data, as it is 
out of the field-of-view of the MOS cameras and just on the edge 
of a pn CCD. 

4.3.3. HD 160529 

The XMM-Newton observation appears more densely populated 
than in the cases of P Cyg and AG Car, but the closest X-ray 
source appears at 56" of HD 160529 and, again, the LBV clearly 
remains undetected. 

4.3.4. qF 362 

The deep Galactic Center Survey data do not show any X-ray 
source in the vicinity of qF 362. The closest X-ray source, 
CXOUGC J174616.6-284909, hes at more than 18", which is 
clearly too far away from the LBV considering the narrow PSF 
of Chandra. 

4.3.5. GAL 024.73+00.69 

A strong flare affected the XMM-Newton observation, reducing 
the useful exposure to only 3.5 ks. This indeed limits the sensi- 
tivity of the observations: only 7 X-ray sources were detected in 
the field-of-view and none appears close of the LBV (they are at 
separations > 7'). 

4.3.6. Pistol Star 

Before the public ation of the full Ga lactic Center survey, 
iMuno et alJ and lOskinoval (|2005b) analyzed the then 

available Chandra exposures of the region, and they report a 
limiting log[Lx/i'BOL] of -8.2 and -8, respectively. 

The deep Galactic Center Survey data do not show any X- 
ray source in the vicinity of Pistol star The closest X-ray source, 
CXOUGC J 1746 15. 2-284957, lies at about 6". The derived limit 
on the log[Lx/i'BOL] is -6, which may see m at odds with pre- 
vious reported results, though the survey of iMuno et al ] d2009h 
uses more datasets. Though no detail on the exact derivation of 
the log[Lx/LBOL] limit is provided in these previous papers, it is 
interesting to note that such limits are found if one uses the ob- 
served flux, not the absorption-corrected one: as the absorption 
in the direction of the Galactic Center is high, this explains the 
two order of magnitude apparent discrepancy in log[Lx/i'BOL]- 
However, if one wants to compare objects scattered all over the 
Galaxy, it is important to correct the luminosities for the inter- 
stellar absorption - indeed, the 'canonical' log[Lx/iBOL] - -7 
of O-stars only applies if this correction is done. Unfortunately, 
this results in a much weaker constraint on the X-ray emission 
of the Pistol star 
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4.3.7. HD 316285 

Data in the vicinity of this cLBV were collected by XMM- 
Newton during 3 observations. For 0112970201 - 0145, it ap- 
pears so far off-axis that it is out of the field-of-view of MOS2. 
In addition, there are a lot of flares in 0205240101 - 0956. This 
probably explains why the limiting count rate derived with only 
the 0112970101 - 0145 dataset is more constraining than the 
limit derived when considering all 3 datasets. We therefore quote 
in Table|2]the former value only. 

4.3.8. GRS 25.5-1-0.2 

This cLBV was observed quite far off-axis (10') in a relatively 
short (7.5 ks) exposure: the sensitivity is thus quite shallow, and 
the derived limit on the count rate is accordingly high. The clos- 
est X-ray source appears at 4' of the cLBV, excluding any con- 
fusion problems. 

4.3.9. WR102ka 

The deep Galactic Center Survey data do not show any X-ray 
source in the vicinity of WR102ka. The closest X-ray sources, 
CXOUGC J174616.6-290115 and J 1746 18. 3-290204, lie at 
about 28". 



4.3.10. LBV 1806-20 

This cLB Vis situated 14" east of the bright pulsar SGR 1806-20. 
While this separation is not a problem for Chandra, the cLB V is 
situated on the PSF wing of this bright pulsar in XMM.-Newton 
instruments. We did a mosaic with all available XMM-Newton 
observations, except the last one (0604090201 - 1785) which 
was taken with a different filter, but the produced sensitivity 
map has no actual meaning. The situation is even worse than 
for W243 as the separation is smaller: the detection is mostly 
hampered by the closeness of the pulsar rather than the available 
effective area. Therefore, since the cLBV seems undetected, we 
estimate the true detection limit by adding a fake source, with a 
similar PSF, at the position of the cLBVPn We changed its am- 
plitude until it was not detectable. We found that an X-ray source 
100 times fainter than the pulsar's emission would be barely dis- 
cernible: this constitutes our flux limit for XMM-Newton, and 
corresponds to an EPIC count rate of 0.01 cts s ' (Table|2]l. 

For Chandra, there is no confusion problem, but the pul- 
sar left a readout trail in the ACIS data. The background was 
thus not measured in an annulus, but rather taken in a circle of 
the same size as the source region and lying on the same read- 
out trail. This resulted in the limit quoted in Table [3] The same 
regions were indeed used for the HRC data, for coherence rea- 
sons. Note that the cLBV is formally detected at Icr significance 
(i.e. not an actual detection with our criteria, see Sect. 3.1) in 
the ACIS data, but this is clearly an artifact due to the readout 
trail. The Chandra observations are better suited to disentangle 
the pulsar from the cLBV, but they are unfortunately not very 
sensitive, thereby not improving the XMM-Newton results. 



Note that the astrometric errors of the 'XMM-Newton datasets are 
very limited as the error on the pulsar's position, found by comparing 
the Simbad coordinates with those in the detection list of the combined 
mosaic amounts to only 0.2". 



4.3.11. Sher25 

Using a single observation of NGC 3603. lMoffat et al.l(l2002h re- 
ports a non-detection of Sher25, without any detail on the limit- 
ing X-ray flux. 

The additional recent observations enable us to enhance the 
sensitivity in the region by an order of magnitude, resulting in 
a tight upper limit on the X-ray emission of Sher25. Note that 
the source is formally detected with Icr significance (i.e. not an 
actual detection with our criteria, see Sect. 3.1) as there are 3.3 + 
3 net counts in the source region. 

4.3.12. [OMN2000] LSI 

Situated in the outskirts of the massive star-forming region 
W51A, LSI shows a formal Icr "detection" (i.e. not an actual 
detection with our criteria, see Sect. 3.1), with an upper limit on 
the observed flux of 1.4 x 10"'^ ergcm"^ s"'. 

4.3.13. GCIRS 16NE, 16C, 16SW, 16NW 

The deep Galactic Center Survey data do not show any X-ray 
source at the position of these objects. The closest X-ray sources 
are GCIRS 33SE and 34SW (see above), as wefl as CXOUGC 
J174540.1-290025,J174540.0-290028 and J174539.7-290029. 
For all sources, we use the same absorbing column to derive the 
limiting flux. 

This region of the survey is the one observed with the largest 
sensitivity, which explains why the limiting photon fluxes are so 
low. However, this is also the most crowded region of the survey: 
the derived limits may thus be somewhat optimistic, although it 
is difficult to quantify by how much. 

4.3.14. 2MASSJ1 7460562-2851 31 9 

The deep Galactic Center Survey data do not show any X-ray 
source in the vicinity of this target. The closest X-ray source, 
CXOUGC J 174604.4-285 139, lies at about 17". 

4.3.15. MCW930 

The position of this cLBV falls in a gap of the pn camera. Only 
MOS data were therefore used to derive an upper limit on its 
X-ray emission. The closest X-ray source appears more than 3' 
away from the cLBV. There is therefore no doubt on the non- 
detection of this object. 

4.3.1 6. 2MASS J1 71 1 0094-39451 74 

Situated in a field containing only a few faint X-ray sources, 
this object remains undetected in the XMM-Newton observation. 
While a flare affected the exposure, especially the pn data, the 
target was observed nearly on-axis, leading to a r ather deep up- 
per lim it. U nfortunately, as for t he other objects of lWachter et al.l 
(1201 Oh and iGvaramadze et al.l (1201 Obi) , no information on the 
other properties of the object are available, preventing the deriva- 
tion of alog[Lx/LBOL]- 

4.3.17. 2MASS J1 841 5965-051 5409 

No bright X-ray source is present in the field-of-view of the 
Chandra observation, and no X-ray source lies closer than 1' 
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away from 2MASS J18415965-0515409. A flux limit of 1.3 x 
10"'^ ergcm"^ s"' was derived. 

4.3.18. 2MASSJ1 6461 734-4508478 

This source is undetected in the XMM-Newton observation. The 
sole X-ray source of the field-of-view appears at 3.6'. As the 
target was observed not far from the on-axis position (off-axis 
angle=2.6'), a good constraint could be put on its count rate. 

4.3.19. 2MASS J 18022233-2238002 

The region encompassing 2MASS J 18022233-2238002 was ob- 
served three times by XMM-Newton. However, two of these 
datasets (0135743001 - 0600 and 0135743101 - 0607) have a 
bug in their attitude files, which could not be corrected (even 
with the help of the XMM-Newton helpdesk). Therefore, no mo- 
saicking was possible for this target. We instead used a reduction 
process similar to the one applied for Chandra data: we first 
derived for each dataset the number of counts appearing in the 
source and background regions by extracting data in adequate 
zones (2 neighbouring circles of 30" radii), then we found the 
actual exposure times for each dataset at the position of the tar- 
get using individual exposure maps, and finally we use the CIAO 
task aprates to get the 90% detection limit on the count rate. 
However, the two additional datasets did not improve the upper 
limit found when using only the 0135742801 - 0600 exposure, 
probably because of the flares affecting the other Rev. 0600 data 
and the fact that the target appears in a gap of the pn camera for 
Rev. 0607. Table|2]therefore quotes the value derived using only 
the first dataset. 

4.3.20. MN53 

This object was observed by both XMM-Newton and Chandra. 
Both observations were centered on the bright IRXS J170849.0- 
4009, and the field contains a few additional faint sources (de- 
tected only with XMM-A^ewfo«). MN53 remains undetected, the 
closest X-ray source appearing at 70" in the XMM-A^ewfo« data. 
Note that the pn data are limited to a small area around the main 
target, so that the flux limit was derived using only MOS data. In 
the Chandra data, MN53 is formally detected at Icr significance 
(i.e. not an actual detection with our criteria, see Sect. 3.1) and 
the flux limit derived on this HRC observation is about 20 times 
larger than that found from the XMM-Newton data of the same 
region. 

4.3.21. MN79 

MN79 appears in a sparsely populated region, containing only 
faint X-ray sources. The closest one in the XMM-Newton data 
lies at 47" of the target and at 15" in the Chandra data, but the 
target itself remains undetected. The two observatories provided 
similar flux hmits. 

4.3.22. MN83 

Situated on the other side of the aimpoint compared to GAL 
026.47-1-00.02, MN83 shows a contrasting behaviour by remain- 
ing totally undetected. The closest X-ray source is at ~ 30", and 
the limiting count rate derived for MN83 is about two orders of 
magnitude below that of GAL 026.47 -1-00.02. 



5. Discussion 

This first LBV survey in the X-ray domain covers half of the 
known L BVs or candidates (or two-thirds of the catalog of 
IClark et a l. 2005b). Despite the heterogeneity of the dataset, the 
derived hmits in observed fluxes are remarkably similar, lying 
between 10"'^ and 10"'"* ergcm"^ s"'. However, as the distance 
and reddening of the surveyed objects vary a lot, so do the con- 
straints on the X-ray luminosity of LBVs: the derived limits span 
more than four orders of magnitude, from ~ 10""' to ~ 7x10"''* erg 
s-i. 

Of all surveyed objects, only four are detected in X-rays, 
corresponding to a detection fraction of one out of four in the 
sample of confirmed LBVs, one out of 9 for the Clark's cata- 
log, or one out of 17 for the full sample of (c)LBVs. Besides 
the two well-known cases of eta Car and CygOB2#12, GAL 
026.47+00.02 and W243 were found to emit some X-rays. Two 
additional cases are formally detected, but the properties of their 
X-ray emission cast doubt on the identification of the high- 
energy source with the (c)LBVs. 

With the results of our survey at hand, what can be concluded 
about the X-ray emission of LBVs? 

First, LBVs, as a class, are clearly not bright X-ray emit- 
ters. Moreover, it seems unlikely that LBVs possess an intrinsic 
X-ray emission of a level comparable to that of O-type stars, 
which are moderate X-ray emitters. While O-stars display a 
tight relation between their bolometric and X-ray luminosities 

(log[Lx/iBOL] 7 with dispersions down to 0.2 dex for clus- 

ters, see clusters studies of [S ana et al. 2006; N aze et al]|201 Ibl 
and discussions in lNazell2009t l2011a). five (c)LBVs m our sam- 
ple have log[Lx/i'BOL] < -8.2, PCyg even showing no emis- 
sion at our sensitivity limit of log[Lx/iBOL] - 9.4! This indi- 
cates a significantly weaker X-ray emission, compared to O- 
stars, for the LBVs as a class. A rather similar conclusion was 
found for Wolf-Rayet stars, especially those of the WC type 
dOskinova et all l2003h . This could have two causes: either the 
local absorption by the wind is too high to let any X-rays gen- 
erated in the deep inner layers escape, like for WC stars, or the 
wind velocity is too small (a few hundred km s~', as is often ob- 
served in LBVs) to produce hot plasma through intrinsic wind 
shocks. While both may well explain the observed lack of X- 
rays in LBVs, it is unfortunately difficult to accurately constrain 
the wind velocity (see e.g. the case of CygOB2#12 discussed 
above) or absorption (see e.g. the case of eta Car, Corcoran 201 1]) 
whatever the phase of LBV activity. Both hypotheses (and a 
combination of the two) therefore remain viable at the present 
time. 

If LBVs are intrinsically X-ray faint, the few cases of X- 
ray detections should have an extrinsic cause. The most obvi- 
ous candidate is binarity. It is indeed well known that wind- 
wind collisions in massive binaries (i.e. O+OB or WR-i-OB) 
can prod uce bright, hard jind variable X-ray emission (for a re- 
view, see Giidel & Naze, 2009). This is certainly the case of the 
Galactic LBV eta Car (see a bove) as well as of the SMC LBV 
HD5980 (iNaze et all l2007l and references therein): the phase- 
locked variations of their bright and hard X-ray emission are a 
definite proof of colliding winds. Such a process has also been 
strongly advocated for Cyg OB2 #12 in view of its spectral prop- 
erties and its luminosity changes (see Sect. 4.1.3 above), though 
a direct evidence for periodicity in the variations is still lacking. 
In addition, the bright and hard X-rays of GAL 026.47+00.02 
are again reminiscent of a wind-wind collision but additional X- 
ray data are needed to unveil the variability properties of the X- 
ray emission and ascertain the nature of this cLBV (Sect. 4.1.4). 
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Finally, W243 does not show any sign of colliding wind X-rays, 
but its X-ray emission may well be that of an O-type companion 
(see Sect. 4.1.2 above). For the four detected (c)LBVs, binarity 
constitutes a plausible explanation, even if other processes, such 
as magnetically-confined winds, cannot be formally ruled out. 

Does that mean that other (c)LBVs are single, as they are not 
X-ray bright? It must be recalled that both Chandra and XMM- 
Newton hav e shown th at X-ray bright CWBs a re the exception , 
not the rule dSana et al .. 2006; Naze, 2009; Na ze et all 1201 Ibh . 
In addition, even in the rare cases of X-ray bright CWBs, the 
X-ray emission may not be as bright and hard thr oughout the 
whole orbit (for a review, see IGiidel&Nazd. 12009). Therefore, 
while a bright, hard and phase-locked X-ray emission may be a 
strong indication of binarity, the contrary is certainly not true. 
The presence of a companion cannot always be inferred from 
X-ray studies, and the single status of the other (c)LBVs of our 
study cannot be definitively established on the sole basis of their 
apparent lack of bright X-rays at the time of the observation(s). 

Our study provides additional information allowing to con- 
strain the presence of an early OB companion as such stars are 
not X-ray dark. They display typical X-ray lum inosities of 10^' 



y ty pi 

to 10^^^ erg s' dHarnden et al.lfl 979; Ku&Chanan 1979l see 



also recent statistical studies of [ Sana et al.i.2006; Antokhin et al.l 
l200i lNazil2009t iNaze etalj|201 Ibh . Even if bright (> lO^^ erg 
s~') wind-wind collisions do not exist in massive systems, the 
intrinsic X-ray emission of the companion may show up, as the 
luminosity ratio is less favorable for the LBV at these high en- 
ergies. Unfortunately, many of our luminosity limits are quite 
high, hence not always able to test this scenario. Only in 5 cases 
(PCyg, AG Car, HD 160529, HD 316825, and Sher 25) are the 
data sensitive enough to probe the range of X-ray luminosities 
of OB stars, taking into account the interstellar absorption of our 
targets. 

If these LBVs possess OB companions, there are basically 
two possibilities: either (1) the companion is quite close to the 
LBV or (2) the companion is relatively distant from the LBV. 
In the first case, the intrinsic emission of the companion may be 
hidden by the strong absorption of the dense wind, but then the 
conditions are favorable to get an X-ray bright wind-wind colli- 
sion. In the second case, no emission from a wind- wind collision 
is expected due to wind dilution, but the intrinsic emission of the 
companion would be easily detectable, as the tenuous (at these 
distances) wind of the LBV cannot hide it anymore. 

To evaluate the limit between these two scenarios, we focus 
on the most constrained case, that of PCyg. The wind density 
p(r) can be written as MjAn r^v(r) where the velocity is assumed 
to follow a beta-law with /3 - 1, i.e. v(r) - Vt»(l - (Rt/r)). In 
this case, the equivalent absorbing column for hydrogen, from a 
position 7?start in the wind can be written as: 



NJr) = ^ r 



p(r)dr — 



-ln(l - 



Using atmosphere models, iNaiarro et al.l ([T997I) and |Nai arrd 
(1200 ih have calculated the stellar parameters of PCyg: R» - 
76 Rq, M ~ 2 X IQ-^ Moyr-i, nHe/nn = 0.30 and = 185 km 
s^^ With these parameters and /?start of 2, 5, 10, 50, and 100 /?», 
we find absorbing columns of 20, 6, 3, 0.6, and 0.3 xlO^^ cm"^, 
respectively. An optically-thin plasma with a temperature of 
0.6 keV and the strength of a late 0-early B star (i.e. 10^' erg 
s^'at a distance of 1.7 kpc), which would suffer both the inter- 
stellar absorption of P Cyg and an absorption by such columns 



" which corresponds to an abundance fraction (in weight) of 
-0.45 



(modelled within Xspec using the vphabs model and a 3 times 
solar abundance for helium) would have been detected as long 
as /?staii > 25/?*. This means that only very close binaries would 
escape detection, and only in the worst cases of orientation (i.e. 
O-star seen through the LBV wind, not through the O-star wind). 
We therefore consider as highly probable that PCyg is not a 
LBV+OB binary. 

In summary, we have information on the distance, reddening 
and X-ray luminosity of 24 objects in our survey (i.e. all studied 
cases except the Wachter et al. and Gvaramadze et al. candi- 
dates). Of these, 4 appear as probable binaries and 5 as probably 
single stars. The status of the remaining objects is unclear: we 
thus have between 4 and 19 binaries and between 5 and 20 single 
stars. Assumi ng a binary fraction of 0.5, as is often suggested for 
massive stars dSana & Evansll20TT I). we modelize the LBV pop- 
ulation result by a random picking of stars using a binomial dis- 
tribution, so that we find the probability to get the different plau- 
sible scenarios (i.e. 4 binaries + 20 single, 5 binaries H- 19 single, 
... 19 binaries + 5 single). The probability to observe our situa- 
tion is found to be 99.91%, for the assumed binary fraction of 
0.5. Only binary fractions below 0.26 or larger than 0.69 would 
lead to a different situation in 10% of the cases. However, such 
very low or very high bina ry fractions seem unl ikely for pop- 
ulations of massive objects dSana & Evansl I20T1I) . Considering 
only the confirmed LBVs yield similar results, though the sam- 
ple is smaller and the statistics therefore less good. In what con- 
cerns the binary fraction, we thus did not find any evidence that 
LBVs are different from other populations of massive objects. 

It must be noted, however, that binarity was proposed to 
play a role in the occurrence of giant eruptions. More specifi- 
cally, these events would be triggered by close interactions be- 
tween the components during some periastron passages. An ex- 
ample would be the 19th-century erupt ions of eta Car, which 
occurred within w eeks of periastron (iKashi & Sokerl 1201 Ot 
'Smith et al.", 2011), and the 20th-century eruptions of HD5980 
(Koenigsberger, 2004; Koeni gsberger et al., 2010). The last case 
of historical giant eruptions is P Cyg, and the impact of a com- 
panion has also been strongly advocated for this star by iKashil 
(2010). While our observations most probably rule out the pres- 
ence of a massive companion, we cannot exclude that a lower- 
mass object lurks in the neighbourhood of PCyg. For iKashl 
(120101) . a late-B companion with 3-6 Mq is sufficient to trig- 
ger massive ejections. Such an object would produce no X-ray 
bright colliding wind emission (since its wind is negligible) nor 
a large amount of intrinsic X-rays: it would thus remain unde- 
tected even in our sensitive observations! The question therefore 
remains open as to whether binarity should be considered as a 
universal trigger of LBV eruptions (see also the discussion of 
alternative triggers in .Smith et al...201 1.) . 

6. Conclusion 

We have performed the first survey of Galactic LBVs and can- 
didate LBVs in the X-ray domain, using a dedicated XMM- 
Newton dataset as well as archival Chandra and XMM-Newton 
data. About half of the known (c)LB Vs (3 1 out of 67) or two- 
third of the (c)LBVs Usted by Clark et al. (2005b, 22 out of 35) 
have now been studied in the X-ray domain. 

The X-ray emission of six objects is detected, but the asso- 
ciation with the (c)LBVs is doubtful for two of them (GCIRS 
33SE and 34W). We are thus left with four (c)LBVs with 
detections in the X-ray range: the well known eta Car and 
Cyg OB2 #12, plus GAL 026.47 -1-00.02 and W243. As it is bright 
and hard, GAL 026.47 -1-00.02 appears rather similar to the two 
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"historical" detections (eta Car and Cyg OB2 #12). These 3 cases 
are well explained, at least quantitatively, by wind-wind colli- 
sions in binary systems. In contrast, the X-ray emission of W243 
appears more modest and it is compatible with the typical emis- 
sion level of OB stars: since the wind velocity of W243 is low, 
the observed emission may be that of an otherwise hidden OB 
companion. 

Excluding these objects, only upper limits were found, with 
values down to Lx = 8 x 10^^ erg s"' and log[Lx/iBOL] - 9.4 for 
PCyg. The constraints on the X-ray emission are in fact particu- 
larly strong in 5 cases (PCyg, AG Car, HD 160529, HD 316825, 
and Sher 25) where even the usual X-ray emission of a poten- 
tial OB companion can apparently be excluded, indicating that 
these stars are probably single and that the intrinsic level of X- 
ray emission is very low in such LBVs. 

Our results (4 detections and 5 strong limits) are fully com- 
patible with binary fractions between 0.26 and 0.69 in LBVs, 
which agrees well with what is found for O and WR stars. To the 
limits of our dataset, the multiplicity properties of LBVs do not 
present significant deviations from other populations of massive 
stars. Binarity may however have an impact on the occurrence of 
giant eruptions, but additional data are needed in order to fully 
test this scenario. 
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